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Abstract 

A review is presented of dilepton and real photon measurements in relativistic 
heavy ion coUisions over a very broad energy range from the low energies of the BE- 
VALAC up to the highest energies available at RHIC. The dileptons cover the invariant 
mass range mi+i- — - 2.5 GeV/c^, i.e. the continmmr at low and intermediate masses 
and the light vector mesons, The review includes also measurements of the 

light vector mesons in elementary reactions. 
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1 Introduction 



Electromagnetic probes - virtual {e~^e~ or ^'^ ^~ pairs) and real photons - are an im- 
portant and powerful tool to diagnose the hot and dense matter produced in relativistic 
heavy-ion collisions (RHI). They are sensitive probes of the two fundamental properties 
that characterize the Quark-Gluon Plasma (QGP) and are predicted by lattice QCD nu- 
merical calculations, the deconfinement of quarks and gluons and the restoration of chiral 
symmetry [T]. 

Their importance has been emphasized time and again since the proposal made by 
Shuryak back in 1978 [2]. Dileptons and photons interact only electromagnetically, and 
therefore their mean free path is large compared to the size of the system formed in nu- 
clear collisions. They are not distorted by final state interactions and once produced can 
escape the interaction region, unaffected, carrying to the detectors information about the 
conditions and properties of the medium at the time of their production. Electromag- 
netic probes are emitted over the entire space-time evolution of the collision, from the 
primordial parton collisions in the early stage till the hadron decays well after freeze-out. 
They are thus very rich in information content but this wealth is also a challenge for the 
measurement and subsequent analysis that need to disentangle the various sources. 

The main topic of interest is the identification of radiation emitted in the form of 
virtual or real photons by the medium in thermal equilibrium. Such radiation is a direct 
fingerprint of the matter formed. Two well distinct sources are considered: 

a) thermal radiation emitted by the strongly interacting QGP in the early phase of 
the collision. The elementary processes involved are the quark annihilation into virtual 
photons [qq — ?■ or the annihilation [qq — )■ 57) and QCD Compton [qg — t- q^y or 

— 57) channels producing real photons. The identification of this signal not only 
serves as a proof of deconfinement but it also provides a direct measurement of one of the 
most basic plasma properties, the temperature, readily given by the inverse slope of its 
expected exponential spectral shape. Theory has singled out dileptons in the mass range 
mi+i- = 1-3 GeV/c^ |2l [3] or real photons in the px range pT = 1-3 GeV/c |4j as the 
most appropriate windows to observe the thermal radiation from the QGP phase. The 
thermal emission rate is a strongly increasing function of the temperature; consequently 
it is most abundantly produced in the early stage when the temperature and the energy 
density have their largest values; the sensitivity to identify it increases the higher the 
initial temperature is relative to the critical temperature of the phase transition. It should 
be easier to observe the QGP thermal radiation at LHC than at RHIC than at SPS. For 
almost two decades this has been a very elusive signal and it is only recently that PHENIX 
at RHIC [5j and NA60 at SPS [6J presented results that could signal the thermal radiation 
from the QGP. 

b) Thermal radiation is also emitted by the high-density hadron gas (HG) in the later 
phase of the collision. The main elementary process here is the pion annihilation into 
dileptons, mediated through vector meson dominance by the p meson (tt+tt" — )• p — )• 
7* —7- l^l^). The HG thermal radiation, controlled by the pole at the p mass of the 
pion electromagnetic form factor, contributes primarily to the low-mass region, [rrii+i- 
< IGeV/c^) around and below the p mass fT]. This component must always be present 
since ultimately the system ends in the hadronic phase and was readily identified in the 
enhancement of low-mass dileptons discovered at the SPS in the mid- nineties jHl |9]. 

The enhancement of low-mass dileptons is one of the main highlights of the CERN 
SPS heavy-ion program. It triggered a huge theoretical activity stimulated mainly by 
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interpretations based on in-medium modifications of the intermediate p meson with a 
possible hnk to chiral symmetry restoration (CSR). (For theoretical reviews see |101 W\\ 
I12j). The SPS results motivated also new experiments aiming at precise spectroscopic 
studies of the vector meson resonances, p, cj and (/>, to explore in-medium modifications of 
their spectral properties (mass and width). This topic is of interest in its own right and 
more so because these modifications could shed light into the behavior of the resonances 
close to the chiral restoration boundary. By definition there are no hadrons in the QGP 
phase. Hadrons are formed in the mixed phase (if there is one) or in the hot and dense 
HG close to the phase boundary; their decay there into lepton pairs can convey valuable 
information about their properties close to the onset of CSR. The p meson with its short 
lifetime, r ~ 1.3 fm/c, and its strong coupling to the vrvr channel is best suited in this 
context, making it the most sensitive signal for in- medium modifications and the best 
probe for CSR. A basic question with no clear answer is how the restoration of chiral 
symmetry is realized. The well established features are that the chiral condensate <qq> 
drops to zero at the phase transition and that in a chirally symmetric world the chiral 
doublets (e.g. p and ai) become degenerate in mass. Unfortunately, this degeneracy 
is practically impossible to establish experimentally. Other implications on the hadron 
spectral properties (mass and width) are uncertain; there is no rigorous connection of 
the chiral condensate to the spectral properties of hadrons. Almost all models predict an 
increase of the width of the p meson at the chiral transition |13| [T^ [151 E] but for the 
behavior of the mass there is no clear picture; in several models the p mass is predicted 
to drop to zero \13\ [T7\ I18j in other models it remains constant |15| . or even increases 
[H]. Under these circumstances experimental insights are very valuable to guide and 
constrain these models. A significant breakthrough has been achieved recently in this 
respect with the precision measurements of NA60 [20j, confirmed also by the upgraded 
CERES experiment [21], which show a significant increase of the width, with no significant 
mass shift, of the p meson at SPS energies. The experimental results seem to be telling 
us that the approach to CSR proceeds through broadening, and eventually melting, of the 
resonances rather than by dropping masses. 

The interest in electromagnetic probes is strengthened by the light vector mesons 
(LVM) studied through their leptonic decays. The special role of the p meson was already 
emphasized in the previous paragraphs; the other vector mesons, uj and (p, can in principle 
reveal the same effects related to CSR as the p. However, the effects are relatively small. 
The ui and (p, with their much longer lifetimes, will decay predominantly outside the 
medium after regaining the vacuum properties; only a small fraction of them will decay 
inside the medium making the effects much harder to observe through direct spectral shape 
analysis. A more promising way to unveil these effects is provided by the simultaneous 
measurement, within the same apparatus, of the (f) meson yield through the e~^e~ and 
K^K^ decay channels. With ~ 2mK even a small in-medium effect on the (j) or kaon 
spectral shapes may have a sizable effect on the (p yield measured through the K^K'' 
decay channel 

The SPS results have also triggered a strong experimental program focused on the 
exploration of in-medium modifications of the LVM in elementary reactions. pA and pho- 
toabsorption reactions at relatively low energies have been mainly studied in this context. 

^The 4> meson has attracted additional attention. The production mechanism in pp vs AA coUisions 
is an open question. Consisting of ss <j) is the hghtest meson with hidden strangeness and can provide 
information about the role of strangeness production in RHI collisions. These are two topics of great 
interest in their own right but outside the scope of the present article. 
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The idea is to avoid the complications due to time variations of the temperature and den- 
sity inherent to RHI collisions by creating a low-momentum particle inside cold nuclear 
matter and observing its decay through dileptons. The disadvantage of course is that the 
predicted effects are much smaller than in high energy nuclear collisions. The difficulties 
in detecting these small effects probably explain the somewhat confusing experimental 
situation. Contrary to the nuclear case where all systems studied show an enhancement 
of dileptons, there is not yet a coherent picture that emerges from the present studies in 
elementary reactions. 

Table 1: Experiments involved in dilepton, photon and light vector meson measurements 
in heavy-ion collisions. The energy is quoted in GeV per nucleon in the lab. system with 
the exception of PHENIX and STAR where the energy is in the c.m.s. 



Experiment 


Probe 


System 


Energy 


Mass range 


Ref 








(GeV) 


(GeV/c2) 




DLS 


e+e~ 


C+C, Ca+Ca 


1.0 


0.05 - 1.0 


[22J 


HADES 




c+c 


1, 2 


- 1 


[23, 24, "25] 


HELIOS-3 


fi+fi- 


p+W, s+w 


200 


0.3 - 4.0 


[9, 26, 






p+Be,Au 


450 




m 


CERES 


e"'"e~, 7 


S+Au 


200 


- 1.4 


t8ll28] 






Pb+Au 


40, 158 




1211 EH |30l BU [32] 


NA38 




p+W, o,s+u 


200 




[331 35^ 36J 


NA50 


fi+fi- 


P+A 


450 


0.3 - 7.0 








Pb+Pb 


158 




[37, '38', '39', 10] 


NA60 




In+In 


158 


0.2 - 5.0 


[6i i20, ^ 42, i3, ^ 


WA80 


7 


S+Au 


200 




m 


WA98 


7 


Pb+Pb 


158 




m 


PHENIX 


e"'"e~, 7 


P+P 


200 


0-4 


laziiiH] 




e+e", 7 


Au+Au 


200 


- 4.5 


laiMlEo] 




(/) K+K- 


Au+Au 


200 




m 


STAR 


p —7- TT+TT^ 


Au+Au 


200 


P 


m 




K+K- 


Cu+Cu, Au+Au 


62.4-200 


4> 


[531 [52 [55] 



The physics potential of electromagnetic probes is confirmed by the wealth of interest- 
ing results obtained so far and by the relatively large number of experiments focusing on 
their study. All the experiments involved in measurements of dileptons, photons and LVM 
in nuclear collisions are listed in Table [l] together with references to their results published 
(or to be published) in the refereed literature. The measurements cover a very broad 
energy range comprising three different energy scales; the DLS at the BEVALAC and 
HADES at the GSI studied dielectron production at low energies of 1-2 AGeV; at CERN 
SPS, low-mass electron pairs were systematically studied by the CERES experiment in S, 
Pb+Au collisions from 40 to 200 AGeV; dimuons over a broader mass range were measured 
by HELIOS-3, NA38/50 and more recently by the high statistics and high mass resolution 
NA60 experiment. There is redundancy in these two energy regimes in the sense that 
at least two different experiments have performed similar measurements with reasonable 
agreement between them giving confidence on the robustness of the results. Given the 
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difficulties of tlie measurements, this redundancy is an important element of the program. 
At the high energies of RHIC, PHENIX and STAR have studied LVM through hadronic 
decays, in particular the (/> — )• K^K^ decay. The study of electromagnetic probes is just 
beginning with first results available from PHENIX on electron pairs. Apart the energy 
scale, the measurements of dileptons are better divided, from the thematic point of view, 
into (i) the low-mass region (LMR, rrii+i- < 1 GeV/c^) where the main contributions are 
Dalitz and two-body resonance decays and where the main topics of interest are CSR and 
in-medium modification of the LVM and (ii) the intermediate mass region (IMR, m;+/- 
= 1-3 GeV/c^) dominated by correlated pairs from semi-leptonic open charm decays and 
Drell-Yan pairs and where the main topic is thermal radiation from the QGP. 

This article is organized as follows. Section [2] discusses the main difficulties of the 
measurement of dileptons and compares the merits of real vs virtual photon measurements. 
Section [3] presents results on p+p and p+A collisions which serve as a crucial reference 
for the nuclear case. Section |4] reviews the results on low-mass dileptons, considering 



first the SPS results from CERES and NA60 in Section 4.1, The RHIC results, mainly 



from PHENIX are presented in Section 4.2 whereas the low-energy results on low-mass 
dileptons at the BEVALAC and GSI are discussed in Section |4.3[ Section [5] presents 
results on the light vector mesons in nuclear collisions. Section [6] contains a discussion on 
IMR dileptons measured at the SPS, first by NA38/50 and more recently by NA60, and 
at RHIC by PHENIX. Section [7| shows results on LVM spectroscopic studies performed 
with elementary collisions. Thermal photons are discussed in Section |8] focusing on the 
attempts of the WA80/98 experiment at CERN and the recent PHENIX results at RHIC. 
A summary and outlook is given in the last Section [9} 



2 Experimental challenge 

The measurement of dileptons is notoriously difficult which explains why the first results 
at SPS and RHIC became available only several years after the beginning of the heavy 
ion program in these machines. There are two main difficulties. The first one is the huge 
combinatorial background of uncorrelated lepton pairs. It arises from the fact that, since 
single leptons do not preserve any information about their parent particle, all leptons 
are paired with all anti-leptons in the same event to form the invariant mass spectrum. 
This background therefore depends quadratically on the particle multiplicity and strongly 
increases as the coverage moves to low-pT leptons. 

In the measurement of e^e~ pairs, the combinatorial background originates mainly 
from vr'^ Dalitz decays and 7 conversions. In an ideal detector with 100% track reconstruc- 
tion efficiency, the combinatorial background would not be a severe problem since the 
tracks from the overwhelming yield of vr" Dalitz decays and 7 conversions form pairs with 
distinctive features (very small opening angle and very small invariant mass). Therefore 
once such a pair is reconstructed, its two tracks can be taken out from further pairing with 
other tracks of the same event thereby considerably reducing the background. However, in 
a real detector, there is a finite track reconstruction efficiency and in particular a minimum 
Pt is required for a track to be reconstructed. Consequently, very often only one of the 
two tracks of a real pair is reconstructed. These single tracks when paired to other tracks 
in the same event give rise to the combinatorial background. 

The combinatorial background (B), can be determined from the yield of like-sign pairs 
or evaluated using an event mixing technique, and then subtracted from the inclusive 
yield of unlike-sign pairs (U) to determine the real signal (S). This subtraction, S = U - B, 
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induces a large statistical and a large systematic uncertainty in S. For example, in the 
measurement of low-mass electron pairs at the SPS, the S/B ratio, before (after) rejection 
has a typical value of ~l/200 (1/20). This implies that the statistical error in S is a 
factor of 20 (6) larger compared to a situation with no background and a 1% systematic 
uncertainty in B results in a factor of 2 (20%) uncertainty in S, AS = 2S (0.2S). 

The combinatorial background is the main factor affecting the quality of dilepton 
measurements and it is therefore imperative to reduce it as much as possible. In the double 
RICH spectrometer of the CERES experiment |56j . this is mainly done by the RICH-l 
detector located before the magnetic field where the original direction of the electrons is 
preserved allowing to recognize conversions and tt^ Dalitz decays by their small opening 
angle. The rejection of close tracks is further supplemented by a doublet of silicon drift 
detectors [31]. After all rejection cuts, CERES reached typical values of S/B = 1/15 -1/20 
in central Pb+Au collisions at 158 AGeV. The NA60 experiment measures dimuons where 
the main combinatorial background originates from charged pion and kaon decays. The 
main rejection tool in NA60 is provided by a precision vertex telescope located close to the 
target. By matching the muon tracks identified after the hadron absorber to the vertex 
telescope the combinatorial background is considerably reduced allowing to reach a S/B 
value of 1/11 in central In+In colhsions at 158 AGeV [20l [57]. The PHENIX detector 
measures electron pairs at mid-rapidity; the present set-up with a strong magnetic field 
starting at the vertex has no rejection resulting in a S/B ~ 1/200 in minimum bias Au+Au 
collisions at ^ys^^ = 200 GeV [59]. The HBD upgrade providing electron identification 
in an almost field free region close to the vertex will considerably improve the situation 
[SHI Eg [60]. 

The second difficulty is the physics background. Photons and dileptons can be emitted 
by a variety of sources and therefore before claiming observation of any new effect, it is 
mandatory to have a thorough understanding of the expected contribution from all known 
sources. Drell-Yan and semi-leptonic decays of charm mesons, produced in the primary 
hard collisions, are the main contributions in the IMR. In the LMR, the physics background 
is dominated by the electromagnetic decays of hadrons (Dalitz decays: vr",//,??' — )• 
Dalitz- vr'^ decay: lo — ?■ l~^l~7r°, and resonance decays: p,ui,(f) — >• /"*"/") which mostly take 
place at a late stage of the collision, after freeze-out. Most experiments have adopted a 
systematic approach, performing precise measurements of the dilepton spectrum in p-|-p 
collisions as the reference for the physics background and p(d)+A collisions to establish 
cold nuclear matter effects. Together they form the basis for identifying any possible 
deviation from the known physics in the hot and dense matter formed in nucleus-nucleus 
collisions. 

The physics potential discussed in the Introduction is in principle as relevant for vir- 
tual as for real photons, since both are expected to carry the same physics information. 
However, their sensitivities are quite different. For real photons the physics background is 
larger than that of dileptons by orders of magnitude, making the measurement of photons 
much less sensitive to a new source [61]. The real advantage of dileptons over photons 
comes from the fact that they are characterized by two parameters, m and pT , whereas 
real photons are fully characterized by one single variable, px ■ 

3 p+p and p+A collisions: the reference measurements 

A precise knowledge of the conventional sources contributing to the dilepton spectrum is 
an essential pre-requisite to unveil effects which are specific to nucleus-nucleus collisions. 



6 



In the LMR, the main dilepton sources that compose the so-called hadronic cocktail, 
are the Dalitz decays, n^, ry, rj' — )• l^l^j, {t:^ contributes only to di-electron measurements) 
the Dalitz-vr*^ decay w — )■ l^l^n^ and the vector meson decays p, cj, (/> — )• /+/^. In di-muon 
experiments there is also a background of correlated pairs coming from hadron decays, 
mainly the p where both pions decay p — ?■ vr+vr" — )• p'^ p~iy^V ^. A weaker contribution 
comes from (p — )• K^K~ — )• p^ p^v^ ^. The contribution of correlated pairs from the 
semi-leptonic decays of open charm mesons, negligible at SPS energies [62] , is an important 
source that has to be included at RHIC energies. Finally hadronic bremsstrahlung that 
contributes to the pair spectrum mainly at very low mass and very low pT was found to 
be significant at very low energies [63] and negligible at SPS energies |64j . 

In the IMR, at SPS energies, the dilepton spectrum is dominated by two sources: 
Drell-Yan pairs and semi-leptonic decays of charmed mesons pairs. Going up to RHIC 
energies, the Drell-Yan contribution becomes negligible whereas the semi-leptonic decays 
of B mesons become significant. 

This reference information is provided by measurements of p+p collisions. The mea- 
surements of p+A collisions allow the study of cold nuclear matter effects. In this section a 
brief review of reference measurements performed at SPS and RHIC is presented together 
with a description of the inclusive dilepton spectrum in terms of these sources. 

3.1 Reference measurements at SPS 

The di-electron sources have been determined at the SPS by the CERES experiment in 
p+Be collisions, a pretty good approximation to p+p collisions. In addition to the inclusive 
e^e~ pair yield, shown in Fig. [!} CERES measured production cross sections and pT 
distributions of vr*^ and r/ via their 77 decay modes, uj via its 7r'^7 decay channel, and fully 
reconstructed Dalitz decays vr'^,r/ — )• 7 e^e~ [27]. These data, complemented with other 
results from NA27 [65], formed the basis of a Monte Carlo event generator, GENESIS, 
used to calculate the inclusive e+e" pair spectrum (see [27j for a detailed description of 
GENESIS and [31j for slight improvements). Mesons that decay into electron pairs are first 
generated according to measured production cross sections, pT and rapidity distributions. 
Mesons are then allowed to decay into e^e~ using known branching ratios [6^ and the 
Kroll-Wada expression |67] with measured form factors for the Dalitz decays |68| |69] or 
the Gounaris-Sakurai expression for the decays of the vector mesons p, a; and (j) into e^e~ 
|70j . The contribution of correlated pairs from the semi-leptonic decays of open charm 
mesons, is estimated using the PYTHIA code ^Tlj. For a meaningful comparison to the 
data, the simulated electron data are finally convoluted with the experimental resolution 
and subject to the same filter as the data in terms of acceptance, pT and opening angle 
cuts. 

The left panel of Fig. [l] shows a comparison of the invariant e^e~ mass spectrum 
measured by CERES in 450 GeV p+Be collisions with the calculated hadronic cocktail [27]. 
One sees that the data are very well reproduced by electron pairs from the known hadronic 
sources. Both the data and the cocktail are normalized to give the pair density per charged 
particle density within the acceptance of the CERES spectrometer. For the cocktail this 
is done by first normalizing to the yield of vr'' which is one of the ingredients of the cocktail 
and then to the charged particle density using the ratio N^^o/Nch = 0.44 [31]. This is a 
convenient normalization for further comparison to other collision systems. If one assumes 
that the relative particle abundances remain unchanged and independent of the collision 
system, then the e^e~ yield is expected to scale with the number of charged particles, and 
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Figure 1: Invariant e^e~ mass spectrum measured by CERES in 450 GeV p+Be (left 
panel) and p+Au (right panel) collisions compared to the hadronic cocktail calculated 
with the GENESIS code. The lines represent the contributions from the known hadron 
decays and the shaded area gives the itlcj systematic error on the summed contributions. 
Statistical (vertical bars) and systematic (brackets) errors are plotted independently of 
each other. The data are not corrected for pair acceptance and the calculated spectrum is 
subject to the same filter as the data, in terms of acceptance, pT and opening angle cuts 

m 

the absolute vertical scale should remain unchanged. This is the case for p+Au collisions 
as demonstrated in the right panel of Fig. [l] that compares the inclusive e^e~ spectrum 
with the same cocktail demonstrating that the cold nuclear matter effects are negligible 
in the dilepton low-mass spectrum at SPS energies. The bands in Fig. [T] reflect the total 
uncertainty in the cocktail calculation which in the mass region m = 0.2-0.8 GeV/c^ is 
dominated by the uncertainties in the branching ratios and transition from factors of the 
7/ and w. These uncertainties are now greatly reduced with the new values derived from 
the low-mass dimuon spectra measured by NA60 in peripheral In+In collisions |43j . 

NA60 performed measurements of p+Be, In and Pb at 400 GeV [72]. The inclusive 
dimuon mass spectrum was analyzed using a similar procedure, using in particular the 
same GENESIS code that was developed by CERES to generate the hadronic cocktail but 
adapted for the case of dimuons. 

The hadronic cocktail for nucleus-nucleus collisions contains the same components 
mentioned above but some of the ingredients need to be adjusted to reflect known changes 
of the relative particle production cross sections, rapidity and pT distributions and their 
dependence with centrality, of the various sources. This is best done by performing inde- 
pendent measurements of the hadrons that contribute to the dilepton spectrum. PHENIX 
does so for the main sources, for the vr'^ — ?■ 77 and similarly for the r] — ?■ 77. For the 
weaker sources where direct experimental information is not available, one can exploit the 
systematic behavior of the data, e.g. using the statistical model that so successfully pre- 
dicts the particle ratios or the systematic behavior of the transverse momentum spectra 
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with mass |31i- NA60 fixes the main sources of the cocktail, for each centrality bin, using 
the measured inclusive dimuon mass distribution itself |20l HT] . 

NA38/50 performed systematic studies of the IMR dimuon yield in 450 GeV p+A 
collisions where A stands for Al, Cu, Ag and W [38j. A very good description of the 
IMR data in all four cases could be achieved by a superposition of Drell-Yan and open 
charm pairs. The event generator PYTHIA was used to calculate the shape of these two 
contributions [7T]. The simultaneous fit of the four data sets, assuming that both cross 
sections scale as A" with a = 1, and imposing the same ratio of the two cross sections for 
the four targets, provided the absolute normalization of the two components with a total 
cc production cross section in p+p collisions at 450 GeV of C7^| = 36.2 zt 9.1/xb. 

3.2 Reference measurements at RHIC 

PHENIX has studied in great detail the dielectron spectrum in p+p collisions at ^/s = 
200 GeV. The inclusive spectrum is shown in Fig. [2] covering the mass range from up to 
8 GeV/c^ [47j. Following the same rational as CERES, a Monte Carlo simulation code, 
EXODUS, was developed for the calculation of the e^e" pair spectrum. The contribut- 
ing sources are the same as those outlined above and the ingredients are similar with 
abundances, pT and rapidity distributions of the various components, taken from other 
PHENIX measurements or from systematics. The main change is that this time the con- 
tribution from heavy flavor (charm and bottom) is substantial. Fitting the data, after 
subtracting the cocktail of the light meson decays, with the charm and bottom spectral 
shapes calculated with PYTHIA, yields a total cross section of = 544 it 39(stat) it 
142(syst) ib200(m-ocieZ) /ib, in reasonable agreement with an independent determination of 
the charm cross section derived from single electron measurements [73]. The comparison 
of the measured e^e~ pair spectrum with the various sources over the entire mass range 
is shown in Fig. [2j 

4 Low-mass continuum in nuclear collisions 

Low-mass dileptons in nuclear collisions have been measured at three different energy 
ranges: 1-2 AGeV at the BEVALAC and GSI, 40-200 AGeV at the CERN SPS and 
y/s~^ = 200 GeV at RHIC. The three experiments that measured lepton pairs at the 
CERN SPS in the early nineties, CERES, HELIOS-3 and NA38, reported almost simulta- 
neously an enhanced production of dileptons (either e'''e^ or ^^^^ ) in S induced reactions 
over a very broad invariant mass range from mu ~ 200 MeV/c^ up to the J /ijj |6lj • Since 
then an enhancement of low-mass dileptons has been observed at all energies and all sys- 
tems studied although it is not clear that the enhancement reflects the same physics in all 
cases. At the two lower energies, measurements were done by two different experiments at 
least, with consistent results. At RHIC energies, PHENIX is presently the only experiment 
capable of measuring low- mass dileptons. 

The excess of low-mass dileptons discovered at the SPS provided a tremendous stimu- 
lus. It triggered a huge theoretical activity mainly motivated by a possible link with GSR. 
It also prompted additional measurements by CERES and motivated a new generation 
experiment NA60 at the SPS. This section reviews the low- mass dilepton results focussing 
first on the SPS results of CERES and NA60, then on the RHIC results of PHENIX and 
finally on the low energy results of DLS and HADES. 



9 




m^e (GeV/c^) 

Figure 2: Inclusive mass spectrum of e^e~ pairs measured by PHENIX in p+p collisions 
at y/s = 200 GeV compared to the cocktail of hadronic sources. The insert is a zoom 
on the low-mass region. Statistical (vertical bars) and systematic (shades) are plotted 
separately. See text for more details [i7] . 

4.1 Low-mass dileptons at the SPS 
4.1.1 CERES results 

The low-mass pair continuum has been systematically studied by the CERES experiment 
at CERN. In addition to the measurements in p-|-Be and p-|-Au at 450 GeV discussed in 
the previous section, CERES measured low-mass electron pairs in S+Au at 200 AGeV [8] 
and Pb-FAu colhsions at 158 AGeV [2ll|29l[3l] and at 40 AGeV [30]. The most prominent 
feature observed in all heavy-ion collisions studied is the enhancement of electron pairs in 
the mass region m = 0.2-0.6 GeV/c^ whereas in p+Be and p+Au collisions the spectrum is 
well described by the known hadron decays. Fig.|3]shows the effect as observed by CERES, 
for the first time in S-|-Au collisions (left panel) and in the last data set in Pb-|-Au collisions 
at 158 AGeV (right panel). The data and cocktail in Pb-|-Au collisions reflect the much 
better mass resolution (Am/m = 3.8%) that resulted from the improved tracking with a 
doublet of silicon drift chambers close to the vertex [73] and the upgrade of the CERES 
detector with the addition of a radial TPC downstream of the double RICH spectrometer 
|75| . With the improved mass resolution a hint of the w and (j) meson peaks is seen for the 
first time in the CERES Pb-|-Au data. The results are presented in the form of electron 
pair production per event and per charged particle within the spectrometer acceptance 
and are compared to the cocktail of known hadronic sources. The yield is clearly enhanced 
with respect to the cocktail. The enhancement factor, defined as the ratio of the measured 
over the calculated yield in the mass range m = 0.2 - 1.1 GeV/c"^, is similar in both cases 
and for the Pb case amounts to 2.45ib0.21(stat)ib0.35(syst)±0.58(decays) (the last error 
represents the uncertainties in the ingredients of the cocktail calculation). CERES has 
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Figure 3: The enhancement of low-mass e~^e~ pairs observed by CERES for the first time 
in S+Au colhsions at 200 AGeV (left panel) [8] and for the last time in Pb+Au collisions 
at 158 AGeV (right panel) [21j compared to the cocktail of known hadronic sources (solid 
black lines), showing also the individual contributions. 

also shown that this excess is mainly due to soft pr pairs and that it increases faster than 
linearly with the event multiplicity |31j . 

This low-mass dilepton enhancement is one of the highlights of the heavy-ion program 
at CERN SPS. Its possible connection to CSR triggered a wealth of theoretical activity 
(for recent reviews see [10^ [TTl I12j). The prime candidate to explain the excess was 
the thermal radiation from the hadronic phase, dominated by the two-pion annihilation 
(vr+vr" — )• p — > e~^e~). This channel contributes a substantial yield of e~^e~ pairs below 
the p mass. However all attempts using the vacuum p properties failed to quantitatively 
reproduce the data. To do so, it was necessary to introduce in-medium modifications of the 
intermediate p meson. Two main venues were used: (i) a decrease of the p-meson mass in 
the dense fireball |76| as a precursor of CSR, following the original Brown- Rho scaling |17j . 
In this scenario, the p- meson mass scales with the quark condensate <qq> and the latter 
drops due to the high baryon density (rather than high temperature) of the medium and 
(ii) a broadening of the p-meson spectral function resulting from the scattering of the p 
meson mainly off the baryons in the dense hadronic medium [IQl \T5[ [T6] . Both approaches 
rely on the high baryon density at mid-rapidity which, at CERN energies, mainly originates 
from baryon stopping, and both achieved good agreement with the CERES data in the 
mass region m = 0.2-0.6 GeV/c2 as illustrated in the left panel of Fig. Q [31]. 

The importance of baryon density prompted the CERES collaboration to measure low- 
mass electron pairs in Pb+Au collisions at the lower energy of 40 AGeV where a higher 
baryon density and a lower temperature are expected at mid-rapidity compared to full 
SPS energy. The results are shown in the right panel of Fig. |4j In line with the theoretical 
expectations, an even stronger, or at least similar, enhancement is observed. The enhance- 
ment factor with respect to the hadronic cocktail, is 5.1ibl.3(stat)ibl.0(syst)ibl.5(decays) 
in the mass range m = 0.2 - 1 GeV/c^. The enhancement is equally well reproduced by 
the two models, dropping mass and broadening of the p, mentioned above [30j. However, 
the errors bars are too large and preclude a more critical assessment. 
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Figure 4: Inclusive e~^e~ mass spectrum measured by CERES in central Pb+Au collisions 
at 158 AGeV (left panel) [31] and 40 AGeV (right panel) [30] compared to the hadronic 
cocktail (thin solid line) . The various components of the cocktail are represented by the 
dotted lines. The figures show also model calculations including the tt+tt" annihilation 
assuming the vacuum p spectral function (dashed), dropping p mass (dash-dotted) and 
in-medium p broadening (solid). 



The success of these two different approaches, one relying on quark degrees of freedom, 
with a direct link to CSR, and the other one based on a many-body hadronic model, 
attracted much debate. The observation that at high temperatures close to Tc, the 
dilepton production rates calculated within the hadronic approach become very similar 
to the qq annihilation rate computed within perturbative QCD raised the interesting 
hypothesis of quark-hadron duality down to relatively low-masses [lUI [77] and provided 
a more direct connection of the widening of the p meson spectral function to CSR. This 
duality will be further discussed below and in the context of the IMR results in Section [6} 

The quality of the CERES data, prior of the upgrade, was insufficient to discriminate 
among the two approaches. A breakthrough was recently achieved primarily with the 
high quality data of NA60 (discussed in the next section) [20] and confirmed by the last 
CERES results with the upgraded spectrometer shown in Fig. [s] [21]. The results are 
compared to dropping mass and broadening calculations [78]. Whereas the two models 
give very similar results for masses m < 0.8 GeV/c^ where the precision of the data is 
still insufficient to discriminate between them, the data between the u and the (j) favor the 
broadening scenario. This is better seen in the right panel, where the excess mass spectrum 
is plotted a-la NA60, after subtracting from the data the hadronic cocktail excluding the 
p meson, and compared to the same model calculations. It is interesting to note that a 
calculation by Kaempfer based on a parametrization of the dilepton yield in terms of qq 
annihilation inspired by quark-hadron duality reproduces the data equally well [79] . 
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Figure 5: CERES invariant mass spectrum of e~^e~ pairs in Pb+Au collisions at 158 AGeV 
(left panel) compared to calculations based on a dropping p mass model (dashed line) and 
on in- medium p spectral function (long-dashed line). The right panel shows the same data 
and the same calculations after subtracting the hadronic cocktail (excluding the p meson) 
ETl. 



4.1.2 NA60 results 

NA60 measured dimuons from threshold up to ~5 GeV/c^ in In+In collisions at 158 AGeV. 
The primary motivation of the experiment was to elucidate the origin of the dilepton excess 
in the LMR discussed in the previous section and in the IMR discussed in Section [6j The 
detector, based on the NA50 spectrometer, was supplemented by a silicon pixel vertex 
tracker |57j that improved the mass resolution to 2.2% at the (j) mass and significantly 
reduced the combinatorial background produced by muons from tt and K decays. The 
vertex tracker also allowed distinguishing prompt dimuons from displaced vertex dimuons, 
a feature essential to clarify the origin of the dimuon excess at intermediate masses (see 



Section 6.1 ). 



The left panel of Fig. [6] exhibits the superb quality of the raw dimuon data in terms 
of resolution and statistics [2(5]. After subtracting the combinatorial background the spec- 
trum contains a total of 360 000 dimuons and the resonances uj and (p are clearly visible 
and resolved; even the ij decay into p'^ p.~ is seen. Whereas the spectrum in peripheral 
collisions is well reproduced by the cocktail of expected hadronic sources -resonance de- 
cays of the r],p,u},(j) and Dalitz decays of the rj,rj ,uj-, this is not so any longer for more 
central collisions. The data show a clear excess of dimuons which increases with centrality 
and is more pronounced at low pair pT ■ The results confirm, and are consistent with, 
the CERES results described previously. The NA60 data quality allows to go one step 
further. The excess mass spectrum, obtained after subtracting from the data the hadronic 
cocktail without the p is shown in the right panel of Fig. |6] (see Refs. [201 HIl EDJ for 
details on the subtraction procedure) . The excess represents the space-time average of the 
p spectral function convoluted with the photon propagator and the Boltzmann factor and 
filtered through the acceptance of the NA60 apparatus. It exhibits a peak at the nominal 
position of the p meson mass, seating on top of a broad structure whose width and yield 
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Figure 6: NA60 data in 158 AGeV In+In collisions. Left panel: invariant mass spectra 
of unlike sign dimuons (upper histogram), combinatorial background (dashed line), fake 
signal (dashed-dotted line) and the resulting signal (lower histogram) [20J. Right panel: 
excess dimuons after subtracting the hadronic cocktail, excluding the p, compared to 
cocktail p (thin solid line), vr+vr" annihilation with an unmodified p (dash-dotted line), 
dropping p mass (dashed line) and in- medium p broadening (thick solid line) |80j . 



increase with centrality. The figure shows also a comparison of the excess with the two 
main models discussed above in the context of the CERES data: in-medium p broadening 
(thick solid line) according to the Rapp-Wambach model |10t [T5t [T6] and dropping p mass 
(dashed line) according to the Brown- Rho scaling |1H I17j. For reference, the figure also 
shows the calculated contributions of the cocktail p , the vacuum p (dash-dotted line), 
and the charm decays (dash-double dotted line). All the curves were obtained after ap- 
plying the NA60 acceptance filter to the theoretical calculations. The comparison clearly 
rules out not only the 7r"'"7r~ annihilation with a vacuum p which was already ruled out 
by the CERES data, but also the dropping p mass scenario. The predictions based on 
the many-body approach of Rapp and Wambach reproduce very well the data for m^+^- 
< 0.9GeV/c^. These early calculations were focussed on the low masses and were not 
intended to reproduce the yield at masses above the p. Recent calculations of the same 
group, that include contributions from multi-pion states, mainly Air states, reproduce the 
entire mass spectrum as shown in Fig. [s] [8T| I82j . 

Other recent calculations reproduce quite well the NA60 excess [831 EH • In Ref. [H3] , 
excellent agreement is found over the entire mass range as demonstrated in Fig. [7j The 
low-mass excess is described also here by the in-medium broadening of the p spectral 
function but at masses above 1 GeV/c^ the thermal radiation from the QGP, and not 
the 47r annihilation, is found to dominate. Ref. |84] provides a very good description of 
the peripheral and semiperipheral mass excess with hadronic rates constrained by chiral 
symmetry arguments and experimental data. 

The excess dimuon mass spectrum, fully corrected for acceptance and reconstruction 
efficiency, with a px cut of 200 MeV/c on the single muon tracks is shown in Fig. ^ The 
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Figure 7: The NA60 dimuon excess in semi-central In+In collisions at 158 AGeV, inte- 
grated over all pT compared to calculations of [83j . The curves represent partial contribu- 
tions and their sum as indicated in the figure. 



figure displays the excess in the LMR discussed here as well as the excess observed in the 



IMR discussed in Section 6.1 The figure also shows the calculations from the three groups 
mentioned above [6l [85j . 
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Figure 8: Absolutely normalized excess dimuon mass spectrum corrected for accep- 
tance and reconstruction efficiency, measured by NA60 in In+In collisions at 158 AGeV 
and compared to theoretical calculations of Renk/Ruppert j83|, Hees/Rapp (SUES] and 
Dusling/Zahed [H]. Both the data and the calculations are subject to a pT cut of 200 
MeV/c on the single muon tracks [Gj I85j. 

The NA60 data pose a new constraint on the theoretical models which have to si- 
multaneously account for the CERES and NA60 results. But the deeper impact of these 
results is their possible relevance to the broader context of chiral symmetry restoration. 
If the system reaches, or is near to, chiral symmetry restoration then the dilepton results 
could be telling us that the approach to such a state proceeds through broadening and 
eventually subsequent melting of the resonances rather than by dropping masses. 
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4.2 Low-mass dileptons at RHIC 

RHIC allows the study of Au+Au collisions up to an energy of ^/s^ = 200 GeV . This 
energy, more than one order of magnitude higher than the SPS top energy of y/s^ = 
17.6 GeV, offers much better conditions for studying the QGP. Higher collision energies 
result in higher initial temperatures or equivalently higher initial energy densities that 
lead to longer lifetimes and larger volumes allowing for a better study of the system before 
it freezes out. Based on the insight gained at the SPS with CERES and NA60, the study 
of low-mass dileptons at RHIC was anticipated to be very interesting. The key factor 
responsible for in-medium modifications of the p meson at SPS energies, as discussed in 
Section |4.1[ is the baryon density or to be more precise the total baryon density. The 
latter is almost as high at RHIC as at SPS. At SPS, most of the baryons at mid-rapidity 
are participant nucleons. At RHIC there is a strong decrease in nuclear stopping but this 
is compensated by a copious production of baryon-antibaryon pairs such that the total 
baryon density is practically the same in both cases [86]. Consequently, the same model 
calculations that successfully reproduce the SPS results predict that the enhancement of 
low-mass electron pairs should persist at RHIC with at least comparable strength [87]. 

PHENIX is the only experiment at RHIC that can measure low-mass electron pairs. 
The mid-rapidity spectrometers have very good electron identification capabilities by com- 
bining a RICH detector with an electromagnetic calorimeter ^88j. The experiment has also 
an excellent mass resolution of ~ 1% at the 4> mass, an essential requirement for precise 
spectroscopy studies of the light vector mesons. However, the experimental set-up is lim- 
ited by a large combinatorial background. The strong magnetic field, starting at the 
collision vertex, causes a limited acceptance of low-momentum tracks in the central-arm 
detectors (which cover each 90° in the azimuthal direction), very often resulting in the 
detection of only one of the two tracks from vr'' Dalitz decays and 7 conversions which 
leads to an overwhelming yield of combinatorial background pairs. With a pT cut of 
200 MeV/c on single electron tracks, the signal to background ratio is of the order of S/B 
~l/200 at rrig+e- = 400 MeV/c^ in minimum bias Au-|-Au colhsions at ^/s^ = 200 GeV 

In spite of that, in a real tour de force, PHENIX mastered the mixed event tech- 
nique for the evaluation of the combinatorial background to an unprecedented precision 
of 0.25%. Fig. [9] shows the first results of PHENIX on e~^e pairs measured in Au-|-Au 
collisions at ^s^^ = 200 GeV after subtracting the combinatorial background HH]. The 
expected yields from the light mesons, calculated using the EXODUS event generator, 
and from the correlated semi-leptonic decays of open charm, calculated using PYTHIA 
and constrained by the measured cc cross section, are also indicated in the figure. At 
very low masses, up to ~100 MeV/c^, the agreement is excellent; at high masses, above 
the (f> meson, it is pretty good also; in particular the continuum in the IMR is fully ac- 
counted by the contribution from charm decays. In the low-mass continuum, from 150 
to 750 MeV/c^ there is a considerable excess of e'^e pairs, with an enhancement fac- 
tor of 3.4ib0.2(stat.)ibl.3(syst.)ib0.7(model) where the last error refiects the systematic 
uncertainty in the total calculated yield. 

The excess is present at all pair px but it is more pronounced at low pair px {pr 



< 0.7 GeV/c) 113. Fig. 10 shows the pT dependence of the excess, after subtracting 
the hadronic cocktail and the charm contribution, plotted for the pair mass range 0.3 < 
rriee < 0.75 GeV/c^ in the form of invariant yield vs mj' — mo {mx is the transverse 



mass niT = \ Pt + "^0)' "^0 is the mean value of dN/dniee in this mass range and the 
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Figure 9: Invariant mass e~^e~ spectrum measured by PHENIX in ^/s~^ = 200 GeV 
minimum bias Au+Au collisions at mid-rapidity. The data are compared to the cocktail 
of expected yields from light mesons and semi-leptonic open charm decays. Statistical 
(bars) and systematic (boxes) errors are plotted separately. The bottom panel shows the 
data to cocktail ratio with the band around 1 representing the systematic uncertainty in 
the cocktail [i^ . 



lower bound of 300 MeV/c^ is chosen to avoid the region of no acceptance due to the pT 
cut, pt > 200 MeV/c, on the single tracks). The distribution exhibits a clear change of 
slope at rriT — rno ~ 0.4 GeV/c and can be well described by the sum of two exponential 
distributions with inverse slope parameters Tl = 92 ±11.4**"* ± 8.4'^J''** MeV and T2 = 
258.3 ±37.3^^*"* ± 9.6*?'** MeV. 

To further characterize this enhancement PHENIX has studied its centrality depen- 
dence. The integrated yield, divided by the number of participating nucleon pairs (Npartf^), 



is shown in Fig. 11 as a function of Npart for two mass windows: < mg+e- < 100 MeV/c^ 
dominated by the tt^ Dalitz decay and 150 < rn-g+g- < 750 MeV/c^ where most of the 
enhancement is observed. The figure also shows the expected integrated yield in the same 
mass windows. Whereas the low-mass yield (bottom panel) is in very good agreement 
with the cocktail, reflecting the expected increase with the pion yield, a very strong cen- 
trality dependence is observed for the low-mass continuum (top panel). The enhancement 
appears concentrated in the two most central bins only, 0-10% where the enhancement 
reaches a factor of almost 8, and 10-20%. 

The PHENIX results in the LMR are intriguing. They appear different from those 
observed at the SPS. The excess is more spread towards lower masses compared to the 
SPS where the excess is closer to the p meson (see right panels of Figs. [5]and|6]). The 
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Figure 10: Invariant niT spectrum of the e~^e~ pair excess (after subtracting the cocktail 
and the charm contribution) measured in the mass range 0.3 < ruee < 0.75 GeV/c^ by 
PHENIX in ^ys^^ = 200 GeV minimum bias Au+Au cohisions at mid-rapidity. The 
soHd hne represents the fit to the sum of two exponential functions shown separately by 
the dashed and dotted lines. Statistical (bars) and systematic (boxes) errors are plotted 
separately [39]. 

centrality dependence is also different. Models that successfully reproduce the LMR results 
of CERES and NA60, based on a broadening of the p meson spectral function fail to 
reproduce the PHENIX results |491 l89l [90] . The results are also intriguing in the IMR but 
these will be discussed later in Section [U 

A new qualitative level in the PHENIX e~^e~ data is expected with the recently in- 
stalled Hadron Blind Detector (HBD) [SSlESEnj. The PHENIX detector was designed 
anticipating that the measurement of low-mass e~^e~ pairs would be feasible with an ap- 
propriate upgrade. In particular, provision was made for the installation of an inner coil. 
When operated in the H — configuration, the inner coil counteracts the main field of the 
outer coil creating an almost field-free region close to the vertex, extending out to ~50-60 
cm in the radial direction. The HBD is located in this field-free region where the pair 
opening angle is preserved and its main task is to recognize and reject electron tracks orig- 
inating from vr'^ Dalitz decays and 7 conversions by exploiting their small opening angle. 
The HBD is a windowless Cherenkov detector operated with pure CF4, in a proximity 
focus configuration. The detector consists of a 50 cm long radiator directly coupled to 
a triple GEM detector which has a Csl photocathode evaporated on the top face of the 
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Figure 11: Integrated yield per pair of participating nucleons as function of Npart compared 
to the calculated yield from the cocktail of light mesons and open charm decays, for two 
different mass windows, < mg+g- < 100 MeV/c^ (bottom panel) and 150 < mg+e- < 
750 MeV/c2 (top panel) from PHENIX pj. 

upper-most GEM foil and pad readout at the bottom of the GEM stack. A large sup- 
pression of the combinatorial background, well above a factor of 20 depending on the cuts 
applied, was demonstrated with p+p data taken during the 2009 RHIC run. A similar 
improvement is expected in the 2010 RHIC run which is devoted to Au+Au collisions. 



4.3 Low-mass dileptons at low energies 

Low-mass electron pairs have also been measured at much lower energies. The DLS 
Collaboration measured low-mass e^e~ pairs in 1 AGeV C+C and Ca+Ca collisions at 



the BEVALAC [22]. The Ca results shown in the left panel of Fig. 12 exhibit an en- 
hancement of low-mass pairs in the mass range m = 0.2 - 0.6 GeV/c^, in comparison to 
transport model calculations based on the HSD (Hadron String Dynamics) approach that 
include all the expected sources: hadron decays, proton-neutron (pn) and pion-nucleon 
(ttN) bremsstrahlung and the pion annihilation {'k'^'k~ — )■ p — )• e'^e~) using the vacuum 
or free p spectral function [63]. The enhancement looks qualitatively similar to the one 
observed at the SPS. But contrary to the situation at SPS energies, all attempts to repro- 



duce it failed. The right panel of Fig. 12 compares the data with the same cocktail using 
this time the p meson spectral function modified in the medium. Part of the enhance- 
ment is explained but the calculations are still a factor of 2-3 below the data [63]. Other 
attempts failed also [QU [92l [93] . This intriguing situation led to the denomination "DLS 
puzzle" that persists for almost a decade and motivated a new experiment, HADES, at 
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the GSI. HADES, was designed and built for systematic studies of low-mass dileptons in 
proton (p-|-p,d,A), pion (vr+pjA) and heavy-ion (C-l-C, Ni-|-Ni, Au-|-Au) induced reactions 
at energies up to 8 AGeV. 
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Figure 12: DLS results in Ca-|-Ca collisions at 1 AGeV j22] compared to HSD calculations 
|63| . The left panel uses the free p meson spectral function whereas the right panel shows 
the calculations with the full p spectral function. 



Several interesting developments have occurred over the last couple of years. From 
the experimental side, first results from HADES are now available \23\ [2^ I25j . They 
are in very good agreement with the DLS results and show an enhancement of low-mass 
dileptons in C-|-C collisions at 1 and 2 AGeV with respect to model calculations. The 
experimental results are therefore on a very robust ground. On the theoretical side, new 
transport calculations based on the HSD approach seem to get close to explaining the ex- 
cess. Using an enhanced bremsstrahlung contribution, in line with recent QBE (one boson 
exchange) calculations (94j , Bratkovskaya and Gassing achieved very good agreement with 
the published DLS and HADES data on G-l-G collisions, as illustrated in the left panel of 
Fig. 13 [95j. The same calculations applied to the Ca-|-Ga system are closer to the data 
(compared to the status of the previous HSD calculations shown on the right panel of 
Fig. 12) but they still underestimate the measured yields. 

Recently, Ref. [96] demonstrated again how crucial it is to have a complete and precise 
knowledge of the elementary reactions for a reliable interpretation of the more complex 
nuclear collisions. In particular, they showed that IQMD calculations are able to reproduce 
reasonably well the HADES C-l-C data at 2 AGeV by lowering the unknown cross section 
pn — )• pnu (cf. right panel of Fig. 13) or/and invoking an in-medium mass decrease of the 
LJ according to m^^ = m^{l — 0.13p/po). 

The measurement of elementary reactions is part of the HADES program and the first 
dilepton measurements on p-|-p and d-|-p collisions at 1.25 AGeV were recently performed. 
The preliminary results are very interesting as they show that the 1/2 (pp -|- pn) spectrum is 
in good agreement with the C-l-C spectrum at 1 AGeV, when both spectra are normalized 
to vr^. The C-l-C data can thus be understood as a superposition of NN collisions and 
there is no compelling evidence to invoke any new source of dilepton or any in-medium 
modification of mesons in this light system and at this low energy [25\ \^7\ . This solves the 
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Figure 13: Invariant mass e~^e~ pair spectrum measured by HADES in C+C collisions at 
2 AGeV |23) compared to model calculations. Left panel: new HSD calculations with an 
enhanced bremsstrahlung contribution and including collision broadening and dropping 
mass effects . Right panel: IQMD simulations with a reduced pn — )• pnu cross section 
using vacuum spectral shapes for the p and lo [UB] . 



DLS puzzle in the sense that there is no excess of dileptons at 1-2 AGeV; the problem is 
reduced to an understanding of the elementary reactions. 

5 Light vector mesons in nuclear collisions 

The previous section demonstrated the particular role of the p in the low-mass dilepton 
spectrum and provided evidence for in-medium effects on its spectral function. Similar 
effects are also predicted for the to and (p mesons. Their observation in nuclear collisions 
is however a formidable experimental challenge. First, the predicted effects are smaller 
compared to the p due to the weaker interaction of the lo and (p with the medium. Fur- 
thermore, with their relatively long lifetimes (assuming that they are unmodified in the 
medium, = 46 fm/c and t^i = 23 fm/c), the u and cf) predominantly decay outside the 
medium after regaining the vacuum properties, with only a small fraction decaying inside 
the medium. Since the measurement integrates over the entire collision's history, this will 
result in a small modification of their spectral shape. Most of the yield will exhibit the 
vacuum spectral shape with only a small fraction showing up as a tail in the low-mass side 
(in the case of dropping masses) or both in the low and high mass sides (in the broadening 
scenario). The identification of these small effects sitting on top of a physics background 
and on top of a usually very large combinatorial background is not an easy task. NA60 
with a mass resolution of 2.2% and a S/B ratio of 1/3 at the <f> mass is best suited for 
these studies at the SPS. First hints of modifications of the uj and (f) have recently been 
reported [6l [80]. At RHIC, PHENIX with its excellent mass resolution of the order of 
1% at the (j) pole should be able to perform spectroscopic studies of the uj and (j) mesons 
once the combinatorial background is reduced with the HBD upgrade. In the meantime, 
PHENIX has attempted a direct spectral shape analysis of the (p meson measured through 
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the K~^K~ decay channel in ^/s^ = 200 GeV Au+Au coUisions [51]. The (f) yield fitted 
with a relativistic Breit-Wigner function convoluted with a Gaussian to account for the 
experimental mass resolution, revealed no significant change in the centroid and width 
values of the (p meson from the PDG accepted values. 

The comparison of the (j) meson decay into K~^K~ and e'^e~ pairs provides an al- 
ternative and powerful tool to search for in-medium modifications. Since the (j) mass is 
close to twice the kaon mass, any in-medium modifications of the (j) oi' the kaons may 
induce a large change in the (j) yield derived from the (j) — )• K^K~ decay mode. This type 
of studies has been carried out at the SPS by NA49, NA50, CERES and more recently 
NA60. The (j) meson was studied via the K^K~ decay channel by NA49 [9^ and via the 
/X decay channel by NA50 |40j in 158 AGeV Pb+Pb collisions. Differences in the (j) 
yield by factors of 2 to 4 were found between NA49 and NA50 in the common my range 
covered by the two experiments, leading to what became to be known as the SPS (j) puzzle. 
Recent reanalysis of the NA50 data [99J shows a smaller difference but the disagreement 
still persists, now at a factor of ~2 (see Fig. 14). There is also significant disagreement 



in the inverse slope parameters derived from the two experiments. The differences are 
too large to be explained by the different rapidity coverage or slightly different centrality 
selection of the two experiments. 



NASO ^+ iL 2000 
NA49K*K' 
CERES K*K' 
CERES e*e 



1 1.2 1.4 1.6 1.8 2 2.2 2.4 2.6 2.8 3 

(GeV/c^) 



Figure 14: Transverse mass distributions of (j) mesons measured at the SPS through the 
K+K- (NA49 and CERES), e+e" (CERES) and /i+zx" (NA50) decay channels in central 
Pb+Pb (NA49, NA50) and Pb+Au collisions (CERES) [99]. 

CERES is the first experiment that measured simulataneously the (j) meson through 
both the e^e^ and K^K^ decay channels within the same apparatus |32j . The measure- 
ments were done in central 158 AGeV Pb+Au collisions. Consistent results were found in 
the two channels within the experimental uncertainties. The K^K~ yield was found in 
good agreement with the NA49 results whereas the e^e~ data were incompatible with the 
original NA50 results. However, this is not anymore the case with the reanalyzed /^^^^ 
data of NA50. Within 1-2 a the CERES e+e~ data are compatible with the new NA50 
results or in other words, the precision of the CERES e^e~ data is insufficient to rule out 
the present level of difference between NA49 and NA50. The (j) puzzle in Pb+Pb is there- 
fore not yet fully resolved. There is room for a substantial difference between the hadronic 
and leptonic decays of the ^ meson, that could result from in-medium modifications of the 
(j) meson [HI] or kaon absorption or rescattering that deplete the (j) meson reconstruction 
in the K+K' channel |1D0]. 

The recent results from NA60 on (j) production via muon and kaon decays in 158 AGeV 
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In+In collisions |44| I1U1| provide additional insight but no coherent picture emerges yet 
from all SPS measurements. In NA60, the yields and inverse slopes from the two decay 
channels are in agreement within errors |lUlj , i.e. the discrepancy seen in Pb+Pb between 
the lJ'~^n~ (NA50) and K^K~ (NA49) decay channels is not seen in In+In. On the other 
hand, the multiplicity dependence of the (p meson yield per participant, integrated over full 
phase space, < 4> > /Npart, is different in the In+In muon channel and the Pb+Pb kaon 
channel. As shown in Fig. 15 < > /Npart 



< (p > /Npart increases monotonically with centrality in 
In+In whereas it seems to saturate for Npart > 50 in Pb+Pb. Furthermore, for Npart > 100 
the (/> meson yield per participant is higher in the muon channel of In+In collisions than in 
the kaon channel of Pb+Pb collisions when compared at the same number of participants. 



:0.08 



A 0.07 

-e- 

V 



0.06 
0.05 r 
0.04 
0.03 
0.02 
0.01 



NA60 (|) 

NA49 i|) ^ KK in Pb+Pb 
NA49 ^ KK in p+p 
NA49 <|) ^ KK in C+C and Si+Si 
CERES (|) -> KK in Pb+Au 
CERES lb ^ ee in Pb+Au 



lJ_ 



50 100 150 200 250 300 350 400 

Npart 



Figure 15: Centrality dependence of the (j) meson yields measured through dileptons and 
kaons at the SPS by NA49, CERES and NA60. The data represent integrated yield over 
the full phase space per participant. Brackets show the systematic errors, for Pb+Pb the 
error bars include also the systematic error [44]. 

The story of the (j) meson at RHIC is similarly unsettled. PHENIX has the capability of 
measuring the meson through the e~^e~ and the K^K~ decay channels at mid-rapidity. 
Preliminary results [102j of the rapidity density per pair of participants in y/s^ = 200 GeV 
Au+Au collisions may indicate a possible larger yield in the dilepton channel compared 
to the kaon one. However, within the large statistical and systematic uncertainties of the 
electron data the two channels yield consistent results. A more definite statement will 
have to await for the improvement in data quality expected with the HBD upgrade of the 
PHENIX experiment. 

We conclude this section by mentioning the STAR results on the p meson reconstruc- 
tion via its tt^tt" decay. STAR reported a downward shift of ~40 MeV in the p meson 
peak position in p+p collisions and a stronger shift of ~70 MeV in peripheral Au+Au 
collisions at ^ys^^ = 200 GeV [52]. The observation has been attributed to in- medium 
modification of the p meson in the dilute hadronic matter of the late stages of Au+Au 
collisions and has been addressed within the two scenarios, discussed previously in the con- 
text of the CERES data. While the collision broadening qualitatively reproduces the trend 
but fails to reproduce the magnitude of the shift |103| , the Brown- Rho scaling succeeds in 
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reproducing the magnitude of the downward mass shift |1U4| . 



6 Intermediate mass region 

Dileptons in the intermediate mass region are particularly interesting. As discussed in the 
Introduction, several early calculations have singled out this region as the most appropriate 
window to observe the thermal radiation from the QGP [^[3|. This section reviews dilepton 
measurements in the IMR at SPS and RHIC energies. 

6.1 IMR dimuons at SPS 

Dileptons in the IMR have been measured at SPS energies by HELIOS-3 (p,S+W at 
200 AGeV) [9], NA38 (p+W, S+U at 200 AGeV) and NA50 (p+Al, Cu, Ag, W at 450 
GeV, Pb+Pb at 158 AGeV) [Ml EH ES] and more recently by NA60 in In+In at 158 AGeV 
[H]. All experiments reported an excess with respect to the expected yield from the two 
main contributions in this mass region, Drell-Yan and semi-leptonic charm decay. An 



example of the enhancement seen by NA50 is shown in Fig. 16 The shape of the excess is 
very similar to the open charm contribution and in fact enhancing the latter by a factor 
of 3.6 nicely accounts for the inclusive spectrum. This is the basis for the hypothesis of 
enhanced charm production put forward by NA38/50 to explain the excess |38j . This line 
of thought was pursued in a more formal way in Ref. ^05j suggesting gluons produced in 
soft baryon-baryon collisions as an additional source of cc pairs. 
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Figure 16: Di-muon invariant mass spectrum measured by NA50 in central Pb+Pb col- 
lisions at 158 AGeV together with the total expected yield (lower solid line) obtained 
from the sum of the Drell-Yan (dashed line) and open charm (dotted line) contributions. 
Adding the latter multiplied by a factor of 2.6 (dash-dotted line) leads to a total yield 
(upper solid line) that closely follows the data ^38j . 

Other possibilities to explain the IMR excess were also explored, like secondary Drell- 
Yan emission |106j and final-state rescatterings that broaden the niT distribution of 
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charmed mesons enhancing the yield in the phase space covered by the NA50 experiment 
|107j . Most of the attention focussed on the possibihty of explaining the excess as ther- 
mal radiation. Gale and co-workers followed a similar line to the one that successfully 
reproduced the low- mass dileptons discussed in Section 4.1 in terms of thermal radiation 
from the hadronic gas. On top of the physics background of Drell-Yan and open charm 
pairs, they systematically considered all sources of muon pairs from secondary meson in- 
teractions that could contribute to the IMR: tttt, np, ttoj, vrai and KK — )■ II and were 
able to reproduce the mass and pT distributions of HELIOS-3 in S+W [108^ I109j and 
NA50 in Pb+Pb collisions (see Fig. 17, left panel) [llOj . On the other hand Rapp and 




Figure 17: NA50 dimuon data compared to calculations of Ref. jllO| (left panel) and 
Ref. |lllj (right panel). The various contributions (open charm, Drell Yan and thermal) 
and their sum are indicated in the figures. The thermal contribution is calculated using 
hadronic (partonic) rates in the left (right) panel. 



Shuryak [111] exploiting the quark-hadron duality, used the pQCD qq annihilation rates 
to calculate the dilepton yield throughout the entire space-time evolution of the collision 
which they described with a simple fireball model. This thermal radiation, added to the 
background (Drell Yand and open charm) contributions, is also able to explain the NA50 
enhancement in the IMR as shown in Fig. 17 (right panel). Only a small fraction of this 
radiation is emitted at the early stages and is associated with the QGP phase. These cal- 
culations also reproduce the transverse momentum dependence of the muon pairs in the 
IMR. Similar results using a similar approach were obtained in Ref. [112] . This ambiguity 
between the hadronic and partonic description of the IMR excess might be resolved with 
the new NA60 results that will be discussed now. 

NA60 measured the dimuon spectrum in 158 AGeV In+In collisions up to 5 GeV/c^ 
and finds also an excess in the IMR |6j. The excess has the same shape as the open charm 
contribution calculated with PYTHIA, the same observation already made by NA50 in 
Pb+Pb collisions and mentioned previously. Using the vertex telescope information, NA60 
is able to go one step further. With a vertex resolution of the order of 10-15^ NA60 can 
distinguish dimuons from a prompt source originating at the vertex and dimuons from 
the semi-leptonic decays of D mesons originating from a displaced vertex. This allowed 
NA60 to establish that the IMR excess is from a prompt source, thereby ruling out charm 
enhancement as explanation of the excess. The yield of this excess increases with centrality 
with a stronger than linear, but weaker than quadratic, dependence on Np^rt- 
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Figure 18: Transverse mass distributions of the excess of muon pairs measured by NA60 
in 158 AGeV In+In collisions for four mass windows in the LMR and the (j) (left panel) 
|4H [HU] and for three adjacent mass windows from top to bottom, 1.16-1.4, 1.4-2.0 and 
2.0-2.56 GeV/c^ in the IMR (right panel) [6]. 



The centrality-integrated rriT distributions in the IMR combined with those from the 
LMR turned out to be very illuminating. The results are shown in Fig. 18 Ignoring 



the steepening of the distributions at low niT for the low mass windows, whose origin 
is not understood, all distributions exhibit a pure exponential shape and can be fit by 
exp{—mT /T). Whereas in the LMR, the inverse slope T monotonically rises with mass, 
reaching values of ~250 MeV at m = 1 GeV/c^, as expected from a hadronic source 
subjected to radial flow, the inverse slope in the IMR is practically independent of mass 
with a relatively low value of ~190 MeV. The inverse slope parameters derived from the 



NA60 data vs. mass are displayed in Fig. 19 
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Figure 19: Inverse slope parameters vs. mass derived from rriT distributions of the dimuon 
excess and hadrons (rj, p,uj,(j)) measured by NA60 in 158 AGeV In+In collisions. The 
bands show the inverse slope parameters for Drell-Yan and open charm calculated with 
PYTHIA 161. 
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To reconcile these two apparently contradictory findings (T increasing with mass in the 
LMR and T independent of mass in the IMR) together with the sharp drop of T in the IMR, 
NA60 concludes that the excess of dimuons in the IMR is the thermal radiation emitted 
early in the collision, before radial flow develops and thus is most likely of partonic origin 
thereby providing an experimental argument that breaks the quark-hadron duality. This 
interpretation is in qualitative agreement with recent calculations using an improved space- 
time description of the fireball that report a substantial or even dominant contribution of 
thermal radiation from the QGP in the IMR [811 \8^ 1113] . It is also consistent with recent 
NA60 measurements showing that the excess dimuons are not polarized as expected for 
the thermal radiation from a randomized source |42J. 



6.2 IMR at RHIC 

The dilepton yield in the IMR looks different at RHIC [12] • Within the experimental 
uncertainties, the IMR yield scales as expected with the number of binary collisions as 



illustrated in Fig. 20 But contrary to the SPS where the IMR yield shows a clear excess 
beyond the expected yield from Drell-Yan and open charm, that is now attributed to ther- 
mal radiation of partonic origin, at RHIC the inclusive yield seems to be exhausted by 
the open charm contribution as estimated from PYTHIA (see Fig. [9]) . There is only room 
for an additional source within the experimental uncertainties. PHENIX also points out 
that the charm contribution calculated with PYTHIA represents an upper limit, the yield 
could be smaller due to in-medium suppression of high pT charm quarks [114] that smear 
the correlation between the correlated cc pair. The lines labeled cc — >• e^e~ (random 
correlation) in Figs, [o] and 20 represent an estimate of the minimum open charm contri- 



bution under the extreme assumption of complete randomness between the cc pairs that 
leaves some more room for additional sources. This is an interesting hypothesis that can 
be addressed experimentally once PHENIX is able to disentangle prompt from off-vertex 
electron pairs with the vertex detector under development [115| . 




Figure 20: Dielectron yield in the IMR scaled with the number of binary collisions Ncoii 
vs. Npart measured by PHENIX in Au+Au colhsions at ^/s^ = 200 GeV [39]. 
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7 Light vector meson spectroscopy in elementary reactions 



The discovery at the SPS of in-medium modifications of the p meson in nuclear cohisions 
and their possible link to CSR, motivated a vigorous program of precision spectroscopic 
studies of the light vector mesons (LVM), p, cj and (p in elementary collisions, proton 
induced and photoproduction reactions |116j . The most prominent results are reviewed in 
this section. 

A list (probably incomplete) of recent and current experiments is given in Table [2} 
It includes measurements of LVM through their e'^e~ decay channel which allows clean 
spectroscopic studies not distorted by final state interactions, but also a few measurements 
where the LVM are identified through hadronic decay modes. The latter benefit from 
considerably larger branching ratios but could be affected by rescattering effects of the 
hadrons in the nuclear medium. Elementary reactions are simpler systems to study. They 
have well controlled conditions as the nucleus is always at T=0 and at constant baryon 
density po as opposed to nucleus-nucleus collisions where the temperature is different from 
zero and where both the temperature and the baryon density vary during the collision. 
The predicted modifications at normal nuclear matter density po are expected to be much 
smaller than in nuclear collisions but they are still large enough to be observable. For 
example, models based on the Brown-Rho scaling [T7j and QCD sum rules [18] predict a 
decrease of the vector meson mass of ~15% at p = pQ. The rho meson is advantageous 
because with its short lifetime it has a large probability to decay inside the nucleus. The 
disadvantage is that one has to disentangle the p and oj yields and possible p—to interference 
effects. The (p meson offers the opposite properties, a narrow and isolated resonance with a 
low probability to decay inside the nucleus. For the lo and (p^ the observation of in-medium 
effects requires measurements down to low momentum values. 

Table 2: Elementary reactions used to study light vector mesons modifications in normal 
nuclear matter. 



Experiment 


Reactions 


Probe 


Ref. 


KEK-PS E325 


p+A— ;> p, w+X 
p+A^ <j)+X 


p,uj ^ e^e~ 
(j) — )• e+e" 


lirnfTTB] 

1119j 


CLAS 


7+A ^VM+A* 


p, (/) — )• e+e^ 


11251 11221 


CBELS A/TAPS 


7+A w+A* 


UJ —7- 7r'^7 


[120, 121, 124J 


TAGX 


7 + yl ^ p + A* 


p —7- TT+TT^ 


jl29til30ul31j 



The experimental results reported so far are very interesting but still insufficient for a 
consistent picture to emerge. The first positive result came from the KEK-PS E325 exper- 
iment that measured the LVM in 12 GeV p+C and p+Cu collisions with an excellent mass 
resolution, ~1% at the (p mass. The measurements were performed in the target rapidity 
region to enhance the decay probability of the virtual photon inside nuclear matter. The 
experiment reported a significant excess of e^e~ in the low-mass side of the uj [1171 1118j 
and (p |119j mesons. The excess is clearly visible in Fig. [2l] (left panel) [118j which shows 
the raw mass spectrum obtained with the Cu target. The combinatorial background shape 
was determined using an event mixing technique. However, since the like-sign spectrum 
was not recorded, the quality of the mixed event technique could not be demonstrated 
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and the background could not be absolutely normalized and subtracted from the raw data 
as usually done. Instead, the normalization was obtained by fitting the mixed event spec- 
trum, together with the relevant hadron decays, to the raw data (see lines in the left panel 



of Fig. 21) 
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Figure 21: Invariant mass spectra, measured in 12 GeV p+Cu collisions by the E325 
experiment at KEK-PS [lIB] . Left panel: raw spectrum with a best fit including the 
individual contributions from hadron decays and the combinatorial background. Right 
panel: spectrum after subtracting the background and the contributions from r] — )• e~^e~ 
7 and w — )• e^e~ vr*^ together with results of a model assuming that the p and u masses 
decrease by 9.2% with respect to their vacuum values at normal nuclear density. 



The right panel of Fig. 21 shows the same data after subtracting the contributions from 
r/ — )• e''~e~ 7 and w — )• e^e~ tt^ and the combinatorial background. The resulting spectrum 
is very well reproduced by a simple model assuming a linear decrease of the vector meson 
masses my with the nuclear density p according to mv{p)/mv{0) = 1 - k^p/po), with no 
in-medium broadening. The best fit to the data, indicated by the solid line in the figure, 
yields a value of k = 9.2% meaning that both the p and oj masses drop by 9.2% at normal 
nuclear density. 

More recently, the same experiment reported the observation of a similar excess in the 
low-mass side of the (j) nieson |119]. The excess was only observed in the lowest pT bin 
covered by the measurement {Pj <1.25). The data required k = 3.4% and also an increase 
of the (f) width by almost a factor of 4. The latter was necessary to generate enough decay 
inside the target nucleus. 

A similar effect was originally reported by the CBELSA/TAPS experiment in the pho- 
toproduction of uj mesons identified through the ir^j decay channel, on Nb and LII2 targets 
|120j . At low w-momenta, a clear excess was shown in the low-mass side of the w meson for 
the Nb case whereas no excess was observed with the LII2 target or at high momenta. The 
data analysis was consistent with a dropping mass of 13%, in agreement with the expected 
predictions of Brown and Rho [17J and Hatsuda and Lee [18]. Subsequent transparency 
ratio measurements on C, Ca, Nb and Pb targets yielded, through comparison to model 
predictions, a strong broadening of the uj meson of 130-150 MeV/c^ HSJ . This result 
together with criticism expressed in the literature [122[ I123j on the way the combinatorial 
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Figure 22: Invariant mass e'^e~ spectrum after background subtraction measured by the 
CLAS experiment in the LVM region together with a fit including the resonance decays 
p, cj, — e+e" and the oj Dahtz decay, using the vacuum mass values [125j . 



background was handled in the data of Ref. |120j , prompted a reanalysis of the same data 
|124j . This reanalysis does not confirm the earlier claim of mass shift. On the contrary, 
the comparison of the data to BUU simulations near the u production threshold (E^ = 
900-1400 MeV) favors the scenario of oj broadening without mass shift. New measurements 
of much higher statistics have been performed to confirm this new finding. 

The CLAS experiment reports no effect in the search for in-medium modifications 
of vector mesons in photo-induced reactions on various targets ^H, C, Fe and Ti over an 
energy range of 0.6 - 3.8 GeV [125 J. The vector mesons were identified through their decay 
into e^e~, with pair momenta of 0.8 - 2 GeV/c, similar to the momentum range covered by 
the KEK experiment. The spectra for C and Ti targets, after background subtraction, are 
shown in Fig. 22 They look very similar to those measured by the KEK group (compare 
to Fig. 21). However, in contrast with the KEK experiment, the CLAS data are very well 
reproduced by calculations based on a transport model using the vacuum mass values of 
the p, cj and (j) mesons |126| . The data analysis yields a very small mass shift of the p 
meson mass, k = 2 it 2%, consistent with zero, and some collisional broadening of the 
p meson width (F = 218 ±15 MeV for the Fe-Ti target, compared to the natural width 
of 150.7 lb 1.1 MeV). The CLAS results are clearly inconsistent with the KEK results. 
They rule out the dropping mass scenarios of Refs. |171 [18] and are consistent with other 
models that predict a broadening of the spectral shape without or with very small mass 
shift [T01 IT271IT28] . 

The TAGX Collaboration measured the p meson, identified via its tt+tt" decay channel, 
in photoproduction reactions on ^H, ^He and at = 0.6 - 1.12 GeV [l29l ESOl 
I131| . In this subthreshold energy region, the p production occurs with the help of Fermi 
momentum. The p's are produced with very low momentum such that they have a large 
probability to decay even inside a small nucleus like ^He. The last analysis of the TAGX 
data shows a mass decrease of the p in ^He, of the order of 45-65 MeV/c^ |131j . in contrast 
to the CLAS results previously mentioned. On the other hand, the ^^C data show mainly 
a broadening of the p, with no mass shift, in reasonable agreement with the many body 
effective Lagrangian approach of Rapp and Wambach [16j or the similar picture of Ref. 

[132]. 
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8 Thermal photons 



Thermal photons are treated in depth in a dedicated article in this volume |133j . Here, 
only a brief account is given since real photons are expected to carry the same physics 
information as virtual photons and are therefore linked to the essence of this article. 

The interest in thermal photons was articulated in some detail in the Introduction. 
The unambiguous identification of thermal photons from the QGP is a very strong signal 
of deconfinement and provides a direct measurement of the plasma temperature. The 
absolute yields, obtained by integrating the emission rate over the space-time evolution of 
the collision, have been calculated by several authors and considerable progress has been 
achieved [H I134j . The calculations point to rather well established features. At tempera- 
tures close to the phase transition, thermal photons from partonic and hadronic processes 
have similar production rates in close analogy to the behavior of dilepton production rates, 
perhaps another manifestation of quark-hadron duality. In an elaborate calculation, the 
transverse momentum range pT = 1-3 GeV/c appears as the most promising window 
where the QGP radiation could shine over other contributions in central Au+Au collisions 
at ^/s^ = 200 GeV [4j. 

The thermal photon yield is relatively small compared to the inclusive photon yield. 
As discussed in Section [2| the physics background for real photons (mainly from tt^ and 
r] decays) is larger by orders of magnitude compared to dileptons (at m > 200 MeV/c^), 
making the measurement of real photons much less sensitive to a new source than dileptons. 
Indeed there is no conclusive evidence for QGP radiation from real photon measurement 
at the SPS. The initial measurements of WA80 [35] and CERES |28j with a S beam at 
200 AGeV yielded only upper limits of ~15% of the total photon yield for any source 
beyond the known hadron decays. In the two experiments the sensitivity was limited 
not by the statistics but rather by the systematic errors, too large to identify a thermal 
photon source which could only be of the order of a few percent of the inclusive photon 
yield. In later measurements, the WA98 experiment observed a photon excess in Pb+Pb 
collisions at 158 AGeV [l6]. The excess is small, only a l-2a effect, and was observed 
only in central collisions (in peripheral events the yield is consistent with the expectations 
from hadron decays) . The WA98 results attracted considerable theoretical interest [H [121 
I135| I136| 11371 11381 1139| . Fig. 23 shows the pT distribution of the excess photons from 
WA98 together with a state-of-the-art calculation that provides a rather good description 
of the data [Ij. The calculation uses a simple fireball evolution of the collisions with a 
reasonable initial temperature of Tj = 205 MeV and includes two main components: (i) 
a contribution of direct photons from primordial hard scattering during the first instants 
of the collision. This contribution saturates the yield at high px when supplemented 
with a transverse momentum k-r broadening (Cronin effect) and (ii) thermal photons that 
dominate the low px part of the spectrum, predominantly emitted from the hadron gas 
{npai gas). It is interesting to note that the latter was calculated using the same many 
body approach that successfully explained the excess of low-mass dileptons discussed in 
Section 4.1 thus providing a consistent description of both the excess of low-mass dileptons 
and the excess of real photons in terms of thermal radiation from the hadronic phase. The 
QGP contribution was found to be relatively small in the entire pT range. 

An interesting and original analysis has recently been presented by PHENIX at RHIC 
[5]. Capitalizing on the idea that every source of real photons should also emit vir- 
tual photons |140j and in order to avoid the huge physics background inherent to real 
photon measurements, PHENIX analyzed e~^e~ pairs emitted in Au+Au collisions at 
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Figure 23: Photon excess in central Pb+Pb collisions at 158 AGeV from WA98 |46j 
compared to a calculation including prompt (pQCD) and thermal (hadron gas and QGP 
emission) photons [1]. 



^Sj^j^ = 200 GeV with low masses (wg+e- < 300 MeV/c^) and high transverse momentum 
{1< Pt < 5 GeV/c). In this restricted kinematic window, PHENIX observes a significant 
excess of e~^e~ pairs beyond the expected yield from the hadronic cocktail of light mesons 
and open charm decays. The e~^e~ invariant mass excess is transformed into a spectrum 
of real photons under the assumption that the excess is entirely due to internal conversion 
of photons. Fig. 24 shows the resulting photon spectrum. The figure includes also the 
direct photon data measured by PHENIX at high pT in Au+Au collisions \A8\ [50] . In the 
Pt range of overlap the two data sets, from dileptons and direct photons, are in good 
agreement within the experimental errors. The figure also shows the baseline spectrum 
measured in p+p collisions. The p+p data are well described by NLO pQCD direct photon 
calculations (see solid lines on the p+p data) and by a power-law fit (dashed line on the 
p+p data) |141j . The Au+Au data on the other hand are clearly above the Taa scaled 
p+p fit (see dashed line on the Au+Au data) and require an additional source. The solid 
lines on the Au+Au data represent the power-law fit to the p+p data, scaled with Taa, 
plus an exponential fit. For this exponential component the fit yields an inverse slope 
parameter which is within errors independent of centrality and found to be T = 221 ± 
19 ± 19 MeV in central collisions. Interpreted as thermal radiation, this would represent 
the temperature of the system averaged over the space-time evolution of the collision and 
would correspond, according to hydrodynamical models to an initial temperature of Tj = 
300 (600) MeV for a thermalization time of 0.6 (0.15) fm/c pl2] . 
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Figure 24: pT distributions of direct photons in p+p collisions (invariant cross section) 
and in Au+Au collisions (invariant yield) for several centrality bins from PHENIX [5]. 
The solid symbols are results from an analysis of low-mass high pT dileptons, whereas the 
open symbols are previous PHENIX results of direct photon measurements |481 [50] . See 
text for more details. 

9 Summary and outlook 

In spite of the experimental difficulties, huge progress has been achieved over the past 
fifteen years in the measurement of electromagnetic observables. Like other aspects of the 
relativistic heavy ion program, the field is largely driven by experiment. 

At SPS energies, all heavy-ion experiments have observed an enhancement of dileptons 
compared to the expected known sources and a coherent picture seems to emerge for 
the entire range of dilepton masses. In the LMR, the enhancement is associated with 
the thermal radiation (mainly vr^vr" — t- p — )■ 7* — t- /^/~) emitted by the high density 
hadronic phase, close to the phase boundary, with strong in-medium modifications of the 
intermediate p. It seems that the onset of chiral symmetry restoration proceeds through 
broadening and subsequent melting of the resonances rather than by dropping masses. 
The enhancement in the IMR is not due to charm enhancement as originally suspected. It 
is also associated to thermal radiation but the temperature derived from the inverse slope 
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of the exponential mx distributions strongly suggests a partonic rather than hadronic 
origin. 

At low energies of 1-2 AGeV, the DLS puzzle could be close to a solution with the 
recent HADES results on elementary p+p and p+d collisions indicating no enhancement 
of low-mass e^e~ pairs in C+C collisions but rather a mere superposition of nucleon- 
nucleon collisions. If confirmed, this implies an onset of the low-mass dilepton enhancement 
somewhere in the energy range of 2-40 AGeV, almost the same energy range where the 
QCD critical point is being looked for in existing (RHIC, SPS) and planned (FAIR, NICA) 
facilities. 

At the high energies of RHIC, first results from PHENIX are now available. The 
observed enhancement in the LMR appears different from the one observed at the SPS 
and may require an additional source. The apparent lack of enhancement in the IMR 
is intriguing. Further insight to elucidate these results will require measurements with 
better signal to background ratio and the ability to recognize displaced vertices, both 
expected with upgrade projects that are being implemented in the PHENIX detector. 
Combining results of low-mass dileptons with measurements of real photons, PHENIX 
identified an excess of direct photons (beyond the expectations based on measurements of 
p+p collisions scaled to the nuclear case with Taa) with an exponential spectral shape. 
Interpreted as thermal radiation, this excess reflects an initial temperature of 300 MeV 
for a thermalization time of 0.6 fm/c. 

A variety of precise spectroscopic studies to identify in-medium modifications of the 
LVM in elementary reactions where the conditions are very well defined, is underway. 
Effects of broadening and dropping mass are reported but there is no coherent picture so 
far and in particular not all rcsTilts seem to be reconcilable with the experience gained in 
nuclear collisions at the SPS. More data are clearly needed. 

Spectroscopic studies of the LVM in nuclear collisions are much more challenging, in 
particular the direct spectral shape analysis of the oj and cj). An alternative method, 
the simultaneous measurement of the (f) through its K^K~ and e"'"e~ decay modes, is 
pursued at SPS and RHIC to provide evidence for the expected in-medium modification 
of the (f). At SPS, no effects are seen in In-I-In whereas the available results on Pb-|-Pb 
are inconclusive. Results from PHENIX at RHIC are also inconclusive due to the large 
uncertainties of the dilepton data. Higher quality results should be available soon. 
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